Spectral non-uniform temperature, non-local heat transfer, and the spin Seebeck effect 
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We present a theory of the spin Seebeck effect driven by subthermal non-local phonon heat 
transfer and spectral non-uniform temperature distribution. The theory explains the non-local 
behavior of the effect arising from the fact that phonons that store the energy (thermal) and the 
phonons that transfer it (subthermal) are located in different parts of the spectrum and have very 
different kinetics. This gives rise to a spectral phonon distribution function that deviates from local 
equilibrium along the substrate and is sensitive to boundary conditions. The theory also predicts 
a non-magnon origin of the effect in ferromagnetic metals in agreement with observations in recent 
experiments. Equilibration of the heat flow out of the substrate to the Pt probe and backwards 
leads to a measurable vertical spin-current produced by the spin polarized electrons dragged by the 
local thermal phonons. We predict specific sample length limits and other dependencies that can 
be probed experimentally, and obtain the correct magnitude of the effect. 



Recently discovered spin-dependent thermoelectric ef- 
fects have opened the doors to a new frontier of spin- 
tronics that merge spin, charge, and thermal physics, 
known as spin caloritronics} ^ The key and most puz- 
zling effect among them is the transverse spin Seebeck 
effect (SSE) in which a thermal gradient in a ferromag- 
net/substrate structure gives rise to spin-currents which 
vary along the length of the sample and are detected via 
the inverse spin-Hall voltage.'^' This effect has been ex- 
perimentally observed using diff eren t ferromagnetic ma- 
terials: metals,!^ semiconductors ,^2121 and insulators.'^ The 
magnitude of the SSE is quantified by the transport coef- 
ficient Sxy = rp , where Vy is the measured ISHE volt- 
age, w is the width of Pt probe, and V^jT = (T2 — Ti)/L, 
where L is the length of the sample, see Fig. [ija) . The 
ISHE voltage is given by Vy = ^^^^^^(jsXs)^, where js is 
the spin current, s is its polarization, Oh is the spin-Hall 
angle of the probe (in Pt, 6h is of the order of one per 
cent) and p is its electric resistivity. The effect is non- 
local, i.e. it depends on the position along the sample 
rather than the local temperature gradient. In addition, 
the size of the sample is usually about 1 cm and such 
a long-ranged information about position can be trans- 
ferred only by phonons, propagating along the insulating 
substrate. 

Here we show that the non-local SSE is a consequence 
of the non-local energy transfer due to sub-thermal dif- 
fusive phonons that are sensitive to the boundary condi- 
tions and give ris e to a spectral non-uniform temperature 
along the sample .^^i^. In addition, we demonstrate that 
while in the insulator the SSE is likely determined by the 
phonon-magnon mechanism, in the conducting ferromag- 
net (e.g., NigiFeiJi' and GaMnAs^), the magnon mech- 
anism is not the only one available. In fact, in recent 
measurements in bilayer F-Pt wire devices, the specific 
geometry excludes long-ranged propagation of magnons 
and leaves only phonons as a source of non- locality.Ifi^ In 



addition, the experiments by Jaworski et al. in Ref. |3| 
were performed on a material with Curie temperature 
Tc — 130 K, considerably lower than the Debye temper- 
ature Oj) = 350 K, and showed that Vjshe oc M at the 
vicinity of the Curie point; i.e. the SSE signal vanishes 
with the magnetization M with the same critical behav- 
ior. This latter fact excludes the magnon mechanism for 
this case. 

The theory of this phonon- electron SSE, which does 
not involve magnons, has three key physical mechanisms. 
The first (i) involves the non-local nature of the sig- 
nal driven by subthermal phonons. In recent measure- 
ments of the SSE in insulatorJl^ the temperature dif- 
ference between thermal magnons and phonons assumed 
in the current theorjl^^ has not been observed, sug- 
gesting the necessity of the concept of spectrally non- 
uniform temperature. This concept originates from the 
fact that in most dielectrics, and also some semicon- 
ductors, the energy transfer is highly non-locaP^'^ be- 
cause of the strong dependence of the diffusion coeffi- 
cient of phonons on frequency: D (uj) cx w^**, if the 
dominant scatterers are point-like^^. In the diffusive 
regime of the experiments the energy relaxation length 
is given by lin (w) = \/ D (w) Tin ('^), where the energy 
relaxation rate is (T, uj) oc T^u. While the ther- 
mal phonons, hw ^ ksT, are equilibrated, the subther- 
mal low-frequency phonons can deviate from the local 
equilibrium due to the rapid low-frequency growth of 
inelastic length 1^ (w) = (T) (T/w)^^^ , which leads 
to non-local kinetics. Even the concept of the tem- 
perature itself is well-defined only for phonons of high 
enough frequency. For the 'thermal' part of the spectrum 
huj > ksT, the distribution function has a Planckian 
form Tlx (w) = (^e^/^sT _ j^'j with a local temperature 
T = T{x). As a result, the phonons which store the en- 
ergy and phonons which transfer it are located in different 
parts of the spectrum. As illustrated in Fig. [ijb), this 
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FIG. 1: (a) Scheme of the SSE experiment. The effect in- 
corporates three key physical mechanisms: (i) subthermal 
phonons whose inelastic length, kn, is of the order of the 
sample size, L, and whose elastic length, l^i, is smaller than 
L, drive the non-local heat propagation along the substrate 
which gives rise to a steady state distribution function that 
deviates from local equilibrium; (ii) equilibration of heat flows 
out of the substrate into the Pt probe and backwards estab- 
lishes the temperature in the probe Tpt 7^ T{x); (iii) the 
different phonon distribution functions in the probe and the 
substrate yield a spin-phonon-drag current, js{x) oc 5T±. (b) 
Spectral phase diagram of phonons as a function of sam- 
ple length. The deviation from local thermal equilibrium, 
Sn{ij,x) = n{ij,x) — n'r(x){^), is illustrated (green curve) for 
X = L/3; here Pph{ijj) oc u}'^ is the phonon density of states. 



spectral separation occurs when l^i (T) ^ 1^ (T) <C L, 
where lei (T) = lei {hio = ksT) cx D{fiio = ksT) cx T^^, 
and kn (T) = {huj = ksT) oc T'^-^. Then the sub- 
thermal phonons whose inelastic length is of the order 
of (but elastic length is much shorter) than the sample 
size L drive the non-local heat propagation along the 
substrate, giving rise to a steady state phonon distri- 
bution function that deviates from local equilibrium for 
huj <C ksT and depends on the position along the sub- 



strate. To describe this non-local effect, it is essential 
to formulate the boundary conditions for the equations 
describing the propagation of the diffusive phonons, see 
below. 

The second (ii) mechanism involves the electron- 
phonon drag. Since the probe is a 'dead end', there is a 
full balance between incoming and outgoing heat fluxes 
such that net heat flux is zero. However, the incoming 
and outgoing fluxes have different spectral distributions, 
because of the inelastic processes in Pt which average out 
the spectrum of the incoming flux, and establish a local 
temperature Tpt{x) different from T[x). The spin drag, 
induced by the phonon flux, is sensitive to the spectral 
content of the phonon distribution function. Hence, de- 
spite zero net heat flux, the spin injection is not zero. In 
the stationary situation, the drag voltage induced by the 
phonons is compensated by redistribution of the electron 
density, so that the total electric current is zero (as well 
as electrochemical potential gradient). However, in the 
presence of a spin polarization, there will be a net spin 
current js — j-\ — polarized along magnetization M: 
unlike its charge counterpart, spin drag is not blocked by 
accumulation of the spin density, which is eliminated by 
SO interaction in Pt. The magnitude of js depends on 
the ratio of the thickness of the ferromagnet, dp, and the 
phonon inelastic scattering length there, if^. The opti- 
mal value of dp for observing the phonon drag SSE is of 
the order of l[^ (T). For too thin ferromagnet, dp <C Z^, 
the phonons cannot effectively transfer their momentum 
to electrons to drag them toward the probe. In the op- 
posite limit, dp ^ ifj^, the phonons equilibrate before 
they reach the region near the probe. (An alternative 
mechanism not considered here is the quantum acousto- 
electric pumpinj^^due to the spectrally non- uniform flux 
of phonons .li^ 

The final (iii) mechanism involves the conversion of the 
spin-current to an electric signal via the ISHE. This con- 
version is most optimal if the thickness of the Pt layer 
is of the same order of magnitude as the spin relax- 
ation length in Pt, which is the case in the discussed 
experiments.'^^ As shown in detail in the following sec- 
tions, the resulting theory gives the correct magnitude 
of the signal, predicts a dependence on magnetization 
SssE OC M, and gives specific temperature and size de- 
pendencies that can be tested experimentally. 



I. RESULTS 

A. Subthermal phonon kinetics 

On Fig. [ijb) we show the spectral phase diagram of 
frequency regions contributing differently to the kinetics 
of phonons. There are two characteristic frequencies, w„; 
and LUbai, determining the propagation of phonons: 

hn (w„;) = L, lei {^bal) = L. (1) 
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For non-local transport we require fiLOni = 
kBT{L/lin{T)y^^^ < ksT. In addition, we will 
not be interested in phonons in the ballistic part of 
the spectrum, w < uJbai- This is legitimate as long 
as ujbai ^ w„/, and determines a maximum length of 
the sample, Lmax, given by the point of intersection 
of the curves lin{uj) and lei{iS), as shown Fig. fltb). 
This gives a temperature dependence Lmax T^'/^. 
For lengths larger than Lmax, the non-local effect is 
due to the fraction of phonons propagating ballistically 
and requires a different formalism, which we will not 
discuss here. The other condition that allows to separate 
thermal phonons from those which produce non-local 
effects is Im (T) ^ L. This gives a minimum length 
of the sample oc T~^/^. For length smaller than 

Lmin even thermal phonons are out of equilibrium and 
spectral separation does not hold. The large ratio of 

kn [T) 1 1 el (T) opens the window Lynin < i < ^max, 

which we are interested in. Hence the sample size should 
be in the range indicated on Fig. [ijb). Estimation at 
T = 10 K (when typical phonon energy is 28 K), gives 



ST(x) 



about few cm and Lv, 



on the scale of mm. 



Recall that the typical size of the sample used for the 
SSE experiments is 1cm. With temperature the width 
of the region of applicability of the theory behaves as 
Lmax/Lmin oc T~^/^ and we expect it to be relevant 
up to 50K. In addition, the temperature is assumed to 
be much smaller than the Debye temperature, T 9d, 
which allows us to ignore Umklapp processes. 

With these specific length restrictions we consider next 
the theory of propagation of diffusive phonons along the 
substrate. Owing to the fact that the low-frequency 
phonons do not primarily interact with themselves but 
with equilibrated high-frequency phonons, one may use 
the following kinetic equation that describes propaga- 
tion of phonons in the insulating substrate, valid for 
huj < ksT: 



Sn {uj, x) 



(2) 



Tin (l^) 

where 5n is the deviation from the local equilibrium 

5n{uj,x) ^ n{uj,x) - nT(x){^)- (3) 

To get a closed set of equations, one needs an equation 
for T {x) , which deviates from the linear behavior due to 
the non-locality of the heat transport. This equation is 
obtained from the continuity of the energy density in the 
system, which in stationary situations reads as V-jg — 0. 

Because of the divergence of D (w) at small cj, the 
heat flux jq is transported by the low-energy part of the 
spectrum. ^^ ^^ The heat current density is given by: 

/•OO 

JQ (a^) = ~ / tiujpph{uj)D{uj)dxn{uj,x)duj, (4) 
Jo 

where pph (ui) oc is the phonon density of states 
(summed over all branches). The integral for jqi^x) di- 
verges and has to be cut off at small frequency (the exact 
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FIG. 2: Correction to the linear temperature dependence as 
a function of position, Eq. (12 1. 



value of the cut off does not enter our results since the 
integral for V • jg converges). Using Eq. the energy 
density continuity equation takes the following form: 



uj'^Sn {uj, x) ppfi (w) = 0. 



(5) 



This equation should hold for all x. Thus, one has to 
solve a system of integro-differential equations. For the 
case of a pulse propagation in an infinite media the non- 
local phonon transport has been studied in Rcf. 11 12 22j 
However, we are interested in a stationary solution in the 
presence of the boundaries. 

On the boundary between the substrate and the heater 
there is a jump in the phonon distribution function, be- 
cause of the abrupt change in the properties of materi- 
als. This leads to a finite thermal boundary resistance 
(Kapitza resistance), which manifests itself through the 
jump ATr- at the contact. If the scattering in the vicin- 
ity of the boundary is mostly elastic, the boundary condi- 
tion consists of conservation of spectral heat current den- 
sity across the boundary. It relates the heat flux through 
the boundary to the jump of the phonon distribution 
function across it. At the left end of the sample (which 
is at heat contact with a reservoir at temperature Ti) it 
takes the following form: 



lei (w) dxn{u!,x)\^^Q 



1 



[n(c^,0)-nT, (w)]. (6) 



The boundary resistance Rsd is assumed to be frequency 
independent. If the heat contacts are in thermal equi- 
librium, RbcL can be related to the thermal boundary 

conductivity hsd — aaTk W^Bd^ where Vg is the 
averaged sound velocity. Note, that in the absence of 
the boundary resistance {Rsd = 0), the locally equi- 
librium distribution function n{uj,x) — nxg(x)i'^) with 
Toix) = Ti + (Ta - Ti)x/L satisfies both the kinetic 
equation and the boundary conditions, so that Sn ~ 
and the non-local effect vanishes. One may easily see 
that 6n oc RBd at not too large values of Rsd- 
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With this it is then possible to write down a closed 
equation for T (x). If the phonon temperature as a func- 
tion of position x is known, the distribution function can 
be obtained from Eq. (|2| and reads: 

n{u,x) ^rS'^'^ {ijj,x)+ (7) 
/ n„ (L - a;<) (a;>) nT(^x') {^) dx' /lin (w) , 



where a;< = min(a;,x') , a;> = max(a;,a:') . The 'source 
term n^"^^ (w, a;) comes from the boundary condition ([Hj 
and is equal to 



where 



(w, x) = g^Z^ [?iTi (i:^) (x) + ut^ (oj) n„ (L - x)] , 

(8) 



IIcj (2^) = cosh 



L 



kn (w) 



, sinh 



L 



kn (w) 



(9) 



and 



Z. 



[2g^ cosh {L/kn (lo)) + (l + gl) sinh {L/Un (c^))] • 

(10) 

Above we have introduced the effective boundary ther- 
mal conductance g^j — {hn {1^) /lei (w)) R]j]i- The second 
term in Eq. ^ describes the process of redistribution of 
phonons along the sample due to diffusion and inelastic 
scattering. Substituting Eq. ([7| into Eq. one gets 
an integral equation for T{x), which can be solved nu- 
merically. This procedure is self-consistent: after finding 
T{x), the distribution function is easily calculated from 
Eq. ([7|. To illustrate the result, we assume the following 
ratios of characteristic lengths of a thermal phonon 



L:k„ {T):hi (T) = 1 : 0.1 : 0.005, 



(11) 



and calculate the correction (5T|| (x) to the linear temper- 
ature behavior 



T(x) 



Ti 



T2-T1 
L 



X + STu (x) , 



(12) 



which is shown on Fig. [2j where we have assumed that 
RbcI = 0.1 and Ti < T2. Although the deviation from 
the linear behaviour is small, it ensures the conservation 
of the energy density of the phonons propagating along 
the substrate. Ultimately, the non-equilibrium correction 
6n{uj,x) is responsible for the SSE effect. On Fig. [IJb), 
the frequency dependence of hujpph (i^) Sn {uj,x) is plot- 
ted close to the colder end (for x = 0.3L). On the hotter 
end, 6n has the opposite sign. 



B. Out of plane spin transport 

After finding the non-equilibrium distribution function 
of phonons 6n (w, x), we next concentrate on the heat and 
spin transport in the vertical direction from the substrate 
to the probe across the magnet. The temperature of 
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FIG. 3: The function h{x) oc ST_i_{x) determining the magni- 
tude and spatial profile of the SSE signal Sxy given in Eq. (|20[|. 



phonons in the Pt probe, Tpt (x), is different from the 
T (x). It is determined by the requirement that the heat 
flux created by non-equilibrium non-local phonons from 
the substrate to Pt is compensated by back-flow flux of 
thermal phonons from Pt to the substrate. The resulting 
temperature difference, dT± (x) = T {x) — Tpt (x), can be 
found from the heat balance equation: 



Lopph (^) (w, x) duj = 0, 



(13) 



where 5N (w, x) 



npt is the difference of the dis- 



tribution function of phonons entering and leaving the Pt 
probe, located at x. Here we neglect inelastic scattering 
of phonons while they pass through the ferromagnetic 
layer {dp < if,^), and have assumed the sound velocities 
to be of similar order in the Pt and the substrate. We 
also assume, that the probe is small enough a <C lei{T), 
so that the influence of the counterflow on the phonon 
distribution function in the substrate can be ignored. 
It is useful to present SN (oj, x) in the following form: 

6N (w, x) = [npi^x) (w) - n-Tpt(a:) i^^)] + 5n (cj, x) . (14) 

Then, the the temperature difference 5T±^ {x) can be cal- 
culated from the equation: 

/•oo 

JTl {x) it (x) cx - / z^Sn {zT, x) dz = (x) , (15) 
-'0 

where h (x) is the dimensionless heat flux supplied to the 
probe by the nonequilibrium phonons. The function h (x) 
can be written in a form of h (x) = ^ (uJbai/T)'^ H (x) , 
where H (x) is a slow function of temperature and bound- 
ary resistance Rpd- Here LUi,ai encodes information about 
scattering of phonons on the disorder and the length of 
the sample. Function H{x) is plotted on Fig. [3] for the 
same sample parameters as before. 

With this we can finally estimate the scale of the 
SSE due to conducting electrons, dragged by out-of- 
equilibrium phonons, in more detail. The guiding idea 
about the scale of the effect follows from the derivation 
of the well known Gurevich formulcP^l^^ for the phonon 
drag. This formula gives for thermoelectric coefficient 
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77 = —j/VT the following expression: 77 oc — ^^^^^ ^3 , 

which is valid when q^l ^ 1 (here = kbT/hvs is the 
wavelength of a thermal phonon). For the dirty case 
qxl ^ 1, the particle current density dragged to the 
probe is given bjE3 



fe (x) CX 



PF 



UpSN {UJ, x) W {ujlei/u) Pph (w) du). 



(16) 

We write r^i , Id for electron-impurity scattering time and 
length in the ferromagnet. The role of electron-impurity 
scattering in Eq. ([I6| is twofold. It enhances electron- 
phonon interaction by slowing the motion of electrons 
(making it diffusive). This is taken into account by the 
form of W{x). On the other hand it diminishes the drag 
effect due to the loss of electron momentum by impu- 
rity scattering. The details of the function W {x) depend 
on the character of the scattering of phonons on defects. 
We assume that phonons scatter on impurities vibrating 
with the lattice^S and W {x) = W-uhix). For tempera- 
tures T <C u/lei, we may use the asymptotic behavior 
Wyb (a; ^ 1) « X. Recalling that the charge current will 
be compensated by an unpolarized backflow charge cur- 
rent from the Pt probe, the total spin current is given 
by the polarized current dragged by the phonons into 
the Pt probe. Finally, we rewrite the expression for the 
spin-phonon-drag current injected into the Pt probe as 



jI [x) (X XmT {T/vsf {T/eof A,i (T) J (x, T) 



(17) 



where A^i (T) — {ksT /ep) {kpl^)^ is a dimensionless 
constant, determined by electrons, X 



M 



level of spin polarization and the dragging factor is 



is the 



where S^y = |e| kpp is a material-dependent constant. 
The factor 2.8 takes into consideration that the en- 
ergy of thermal phonon is 2.8kBT. Assuming that in 

we obtain 



"CTTl, 



Pt, p = 0.9 fifl ■ m and kp = 10' 

si^J ~ 30^V/K. Function H (x) is positive at the cold 
end, meaning the dragging force pushes electrons towards 
the magnet there, according to Eq. (19). Note that (i) 



although the electron-phonon drag is proportional to a 
high power of temperature, see Eq. (17), the final result 



for the SSE coefficient is only weakly temperature de- 
pendent, Sxy OC T. It comes out as a result of the strong 
dispersion of the phonon scattering time in the substrate. 
Although function H{x) in Eq. (20) is also temperature 
dependent, this dependence comes only from the non- 
locality of the phonon collision integral in energy, and 
is relatively weak. Another important property of this 
function is that (ii) it's spatial dependence varies with 
temperature rather slowly. This is because the phonons 
which contribute mostly to the non-local effect have in- 
elastic scattering length of the order of the sample size. 
Varying the temperature mainly results in the shift of the 
relevant phonon energy a;„/, so that the corresponding 
length scale lin{^ni) remains the same. The observations 
(i), (ii) stress the importance of the strong dispersion 0/ 
the phonon scattering. 

Note that the sign of the Seebeck follows the sign of 
H{x) and is positive at the cold end. Taking Or = 
0.0037, ep/ke = lO^X, Od = SSOif, kple^ = 10 and 
the ratio of characteristic lengths as in ( [TT] ) , we find the 
magnitude of the effect at \{)K to be 5 ^ l/iV/Kx^A/. 



II. DISCUSSION AND CONCLUSION 



J{x,T) 



z^SN {zT,x)dz. 



(18) 



Since the spectral densities of the energy and the charge 
currents are proportional to different powers of the 
phonon frequency, the electronic drag due to phonons 
is possible even when the net energy flow is zero. The 
contribution to J in Eq. ( 18 ) arising due to the temper- 



ature difference 5Tj_ (x) between the substrate and the 
Pt probe (first term in Eq. (14)) is dominant. In other 
words, while the non-locality of the effect along the sam- 
ple is carried by the low-frequency phonons, the dragging 
force generating the spin-current is produced by the ther- 
mal phonons. As a result of this intricate joint effort by 
the phonons in different parts of the spectrum, one gets 
(restoring units): 

j: (x) = Xm {ks/hf (T) {T/eof {T/v^f 5T^ [x) . 

(19) 

Finally, for the magnitude of the SEE, Sxy = 
-2{\e\p/h)eHjll^xT, and recalling that 5Ti_{x) cx 
—Th{x), we obtain: 

Sxy = ^H^i^'^e/ (T) XMkphi {60/2.8} H (x), (20) 



In this work we have discussed the main ingredients of 
the phonon dynamics in the substrate that allows to un- 
derstand the spatial profile of the SSE signal. As we have 
shown, to explain the non-local effect, i.e. its dependence 
on the position of the probe along the substrate, one 
must consider the spectral non-uniformity of the phonon 
distribution function which can be interpreted as spec- 
trally non-uniform temperature. A key aspect of the 
non-locality is the explicit introduction of the boundaries 
into the equations describing the propagation of diffusing 
phonons. 

In addition, we have presented a scheme of the non- 
magnon mechanism in the case when the ferromagnetic 
element of the device is conductin^^H and obtain the cor- 
rect magnitude of the effect. Furthermore, the spatial 
profile of the SSE signal, presented in Fig. |3j is very sim- 
ilar to the one shown as a 'universal' profile on the Fig. 2f 
of the Ref. '2', Although the phonon kinetics at tempera- 
tures comparable with 9]j is strongly modified by Umk- 
lapp processes, the measured proportionality between the 
SSE signal and the magnitude of the magnetization in 
Ref. m clearly indicates that near Tc ~ 130 K the effect 
is still dominated by the fiux of the spin-polarized elec- 
trons, instead of the magnon-mediated spin torque. We 
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believe that the difference between the data presented in 
Figs 2 and 3 of Ref. [3]- in particular, the difference in the 
behavior near the Tc, - supports this picture. Two dif- 
ferent samples demonstrate drastically different temper- 
ature behavior. The sample which is thicker and grown 
on a substrate of a better quality has larger peak value of 
both Sxy and thermopower a^-x and also much faster de- 
cay of Sxy at approaching Tq- The stronger thermopower 
observed in the thicker sample demonstrates that in this 
sample phonons lose momentum mainly in collisions with 
electrons, while in the thinner sample, their scattering on 
the defects is more efficient. However, due to strong sen- 
sitivity of the phonon distribution function at the F-Pt 
boundary to the ratio dp/lfn (T), in the thicker sample 
the SSE decays with temperature much faster than in 
the thinner one. As we have already discussed, at large 
dp/lfn (^) phonons equilibrate before they reach the 
probe. Indeed, in the thicker sample (more than three 
times thicker than the thinner one) the effect was not 
even resolved near the Curie temperature within the ac- 
curacy of the measurement. This suggests the need to 
study the dependence of the SSE signal on the thickness 
of the magnetic sample in otherwise identical conditions, 
i.e., keeping the properties of the insulating substrate 
and semiconductor/substrate boundary the same. 



Acknowledgements 

The authors are grateful for useful discussion with O. 
Tretiakov, J. Heremans, R. Myers, and G. Jakob. We 
acknowledge the support by grants ONR-000141110780, 
NSF-DMR-1006752, NSF-DMR-1105512, NHARP, and 
the Paul and Tina Gardner fund for Weizmann-TAMU 
collaboration. 



Author contributions 

K. S. T. carried out theory calculations with assistance 
of A. M. F and J. S. All three authors contributed equally 
to the project planning, writing the text and analyzing 
the results. 



Competing financial interests 

The authors declare no competing financial interests. 



* Electronic address: |tikhonov@physics . t amu . edu| 

^ Uchida, K. et al. Observation of the spin seebeck effect. 

Nature 455, 778-781 (2008). 
^ Jaworski, CM. et al. Observation of the spin-seebeck 

effect in a ferromagnetic semiconductor. Nature Mat. 9, 

898-903 (2010). 
^ Jaworski, C. M. et al. Spin-seebeck effect: A phonon driven 

spin distribution. Phys. Rev. Lett. 106, 186601 (2011). 

* Uchida, K. et al. Spin seebeck insulator. Nature Mat. 9, 
894-897 (2010). 

^ Adachi, H. et al. Gigantic enhancement of spin seebeck 
effect by phonon drag. Applied Physics Letters 97, 252506 
(2010). 

® Uchida, K. et al. Long-range spin seebeck effect and acous- 
tic spin pumping. Nature Materials 10, 737-741 (2011). 
Flipse, J., Bakker, F. L., Slachter, A., Dejene, F. K. & 
van Wees, B. J. Direct observation of the spin-dependent 
Peltier effect. Nature nanotechnology 7, 166-8 (2012). 

* Costache, M. V., Bridoux, G., Neumann, 1. & Valenzuela, 
S. O. Magnon-drag thermopile. Nature materials 11, 199- 
202 (2012). 

^ Bauer, G. E. W., Saitoh, E. & van Wees, B. J. Spin 

caloritronics. Nature materials 11, 391-9 (2012). 
^° Valenzuela, S. & Tinkham, M. Direct electronic measure- 
ment of the spin hall effect. Nature 442, 176-179 (2006). 
Levinson, Y. Nonlocal phonon heat transfer. Solid State 
Communications 36, 73 - 75 (1980). 

Levinson, Y. Nonlocal phonon heat conductivity. Sov. 
Phys. JETP 52(4), 704 (1980). 

Agrawal, M. et al. Magnon-phonon coupling unmasked: 
a direct measurement of magnon temperature. arXiv 
preprint arXiv: 1209.3405 (2012). 



Xiao, J., Bauer, G. E. W., Uchida, K., Saitoh, E. & 
Maekawa, S. Theory of magnon-driven spin seebeck effect. 
Phys. Rev. B 81, 214418 (2010). 

Ziman, J. Electrons and phonons: the theory of trans- 
port phenomena in solids (Oxford University Press, USA, 
2001). 

Levinson, Y., Entin-Wohlman, 0.,W61fle , P. Acoustoelec- 
tric Current and Pumping in a Ballistic Quantum Point 
Contact. Phys. Rev. Lett. 85, 634 (2000). 
Tretiakov, O., Tikhonov, K., Sinova, J. & Finkels'tein, 
A. M. unpublished. 

Zhang, S. S.-L. & Zhang, S. Spin convertance at magnetic 
interfaces. Phys. Rev. B 86, 214424 (2012). 
Peierls, R. Zur kinetischen Theorie der Warmeleitung in 
Kristallen. Ann. Phys. (Leipzig) 3, 1055 (1929). 
Pomeranchuk, I. On the thermal conductivity of di- 
electrics. Physical Review 60, 820 (1941). 
Herring, C. The role of low-frequency phonons in thermo- 
electricity and thermal conduction. Proc. Int. Coll (1956). 
Wilson, T. & Schaich, W. A calculation of nonlocal phonon 
heat transfer. Solid State Communications 50, 3-6 (1984). 
Swartz, E. & Pohl, R. Thermal boundary resistance. Re- 
views of Modem Physics 61, 605 (1989). 
Gurevich, L. The thermoelectric properties of conductors. 
Zhurnal Eksperimentalnoi I Teoreticheskoi Fiziki 16, 193 
(1946). 

Blatt, F., Schroeder, P., Foiles, C. & Greig, D. Thermo- 
electric power of metals (Plenum Press New York, 1976). 
Lifshitz, E., Pitaevskii, L. & Landau, L. Physical kinetics, 
vol. 60 (Pergamon press Oxford, 1981). 
Sergeev, A. & Mitin, V. Effect of electronic disorder 
on phonon-drag thermopower. Phys. Rev. B 65, 064301 



7 



(2001). 

Pippard, A. B. Ultrasonic attenuation in metals. Philos. 
Mag 46, 1104 (1955). 

The recent observatiorP^ of the magnetic proximity effect 
in FM/Pt contact suggests the possibility of another chan- 
nel to contribute to the SSE voltage: anomalous Nernst ef- 
fect. However, it is less universal than the electron-phonon 



drag and in any case must rely on the mechanism which we 
propose for generation of the non-local signal: subthermal 
phonons with energy dispersive diffusion, which is quite 
universal. 

S. Y. Huang et al. Intrinsic Spin-Dependent Thermal 
Transport. Phys. Rev. Lett. 107, 216604 (2011). 



